Microscopic derivation of the Jaynes-Cummings model with cavity losses 
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In this paper we provide a microscopic derivation of the master equation for the Jaynes-Cummings 
model with cavity losses. We single out both the differences with the phenomenological master equa- 
tion used in the literature and the approximations under which the phenomenological model correctly 
describes the dynamics of the atom-cavity system. Some examples wherein the phenomenological 
and the microscopic master equations give rise to different predictions are discussed in detail. 
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I. INTRODUCTION 

The Jaynes-Cummings (JC) model is the fundamen- 
tal model for the quantum description of matter-light 
interaction It describes the dynamics of a two-level 
atom strongly interacting with a single mode of the quan- 
tized radiation field in the rotating-wave approximation 
(RWA). This model has been extensively studied in the 
past three decades 0, Purely quantum effects pre- 
dicted by the model, such as Rabi oscillations and col- 
lapses and revivals of the atomic inversion operator, have 
been observed in the experiments both with microcavi- 
ties [3, @] and with trapped ion systems [7] . 

In cavity quantum electrodynamics, the JC model de- 
scribes the strong-coupling regime of micromasers and of 
one-atom lasers. A realistic description of these systems, 
however, must take into account the photon losses due 
to imperfect reflectivity of the cavity mirrors. Cavity 
losses have been described in the literature by means of 
a phenomenological master equation of the form 



-i [Hjcp] + 7 apa 



If If 
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(1) 



where p is the density matrix of the atom-cavity sys- 
tem, Hjc is the Jaynes-Cummings Hamiltonian, and 7 
represents the rate of loss of photons from the cavity, 
and where we have put h = I. The second term on the 
r.h.s. of Eq. H]) has been derived microscopically in the 
framework of a different physical problem, i.e., the one 
wherein the system is given by the quantized cavity mode 
only, losing excitations because of its interaction with the 
surrounding electromagnetic field in the vacuum state Q . 

In this paper we will show that, when the system con- 
sists not only of the cavity mode, but also of the two- 
level atom interacting with the cavity, the fully micro- 
scopic derivation gets more complicated, and in general 
the master equation for the atom-cavity system is not of 
the form given by Eq. ^ . We will show that only under 
certain conditions the phenomenological master equation 
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coincides with the microscopic one. These conditions are 
typically met in the cavity QED experiments, and this ex- 
plains the extensive and successful use of Eq. ([T|) for the 
description of the dynamics of the JC model with losses. 
The knowledge of the limitations of the phenomenologi- 
cal model, however, allows us on the one hand to point 
out some misuse of the model, e.g., for the study of the 
damping of highly excited quasi-classical states f8l|. On 
the other hand it paves the way to the correct description 
of photon losses in the case of structured reservoirs, e.g. 
for photonic bandgap cavities. Moreover, our derivation 
brings to light the microscopic processes which occur in 
the open system dynamics, and therefore how decoher- 
ence and dissipation come into play. 

The paper is structured as follows. In Sec. II we give 
a short review on the JC model and on the usual phe- 
nomenological description of losses in this model. In 
Sec. Ill we present the microscopic derivation of the JC 
model with cavity losses, for the general case of a T tem- 
perature electromagnetic reservoir. In Sec. IV we de- 
scribe the dynamics for the case of one initial excitation 
in the atom-cavity system, we solve the master equation 
and we compare it with the solution of the phenomeno- 
logical model existing in the literature. In Sec. V we 
briefly discuss some situations in which the use of the 
microscopic model could lead to substantial differences 
compared to the phenomenological model, suggesting ex- 
perimental situations in which such differences could be 
observed, and finally we present conclusions. 



II. THE JAYNES-CUMMINGS MODEL AND 
THE PHENOMENOLOGICAL DESCRIPTION OF 
LOSSES 

Since its introduction, in 1963 [1], the JC model has 
been one of the most used models for the description 
of radiation-matter interaction in quantum optics [4], in 
particular in cavity quantum electrodynamics [3| , and in 
ion traps 0] ■ In this section we briefly review some of its 
features. 

Let us consider a two-level atom and denote by \g) and 
|e) the ground and excited state, respectively. The energy 
separation between the two states is given by hujo, with 
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ojQ the Bohr frequency. The resonant interaction between 
the atom and a mode of the electromagnetic field, in the 
RWA and in units of H, is described by the following 
Hamiltonian 



H 



JC 



+ uJo a^a + n (a(T-|_ + aV_) 



(2) 



where (a) is the creation (annihilation) operator of the 
mode, cr_ = \g) (e|, a+ = |e) {g\, and ^ |e)(e| - \g){g\. 

It is straightforward to show that the total number 
of excitations in the atom-cavity system, given by J\f = 
^fltfl + (T2 + 1/2), is a constant of motion. This allows to 
diagonalize easily the Hamiltonian Hjc- One finds the 
following eigenstates and eigenvalues 



\EN,±)^^{\N,g)±\N-l,e)), 



(3) 



for > 1, while the ground state and the corresponding 
energy eigenvalue are 



\Eo)^\0,9) 



— — 



(4) 



respectively, where \N,i) — \N)\i), with i — e,g, indi- 
cates the tensor product of the Fock state |iV) and the 
electronic states \i). 

From the eigenstates and the eigenvalues of Hjc one 
can calculate the evolution of the system given any initial 
conditions. Well known examples of system dynamics 
are Rabi oscillations of the atomic state population, and 
collapses and revivals of the oscillations when the mode 
is initially in a coherent state . 

In cavity quantum electrodynamics the main source 
of dissipation originates from the leakage of cavity pho- 
tons due to imperfect reflectivity of the cavity mirrors. 
A second source of dissipation and decoherence, namely 
spontaneous emission of photons by the atom, is mostly 
suppressed by the presence of the cavity, and therefore 
its effect is usually neglected. 

The dissipative dynamics of the quantized modes of the 
radiation field inside the cavity in absence of the atom, 
i.e., when the atom is not inside the cavity, can be de- 
rived microscopically assuming that the cavity modes are 
coupled with the electromagnetic field outside the cavity, 
which represents a reservoir at temperature T Q. The 
typical approach to the description of the losses in the JC 
model consists in assuming that the presence of the atom 
inside the cavity does not modify strongly the mecha- 
nism of cavity losses, which therefore can be modelled 
by the above mentioned master equation. It is worth 
stressing further, however, that this approach is purely 
phenomenological, since it does not take into account, in 
the microscopic derivation, the presence of the atom in- 
side the cavity. The phenomenological master equation 



has the form 

P = -i[Hjc,p] 

+ 7[n(wo) + l] 
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(5) 



with n{ujQ) the average number of quanta of the reser- 
voir in the mode of frequency luq, and 7 the rate of loss 
of cavity photons. We note that for a zero-T reservoir 
the master equation above reduces to the one given by 
Eq. ((T|). Equations (P) and ([5]) have been assumed to be 
valid in most of the earlier studies dealing with the JC 
model with losses, for example in [1 i, E, El, [S IH, Hi • 

The use of the master equations given by Eqs. ([T]) 
and ([5]) has also been motivated by the fact that in many 
cavity QED experiments the atoms fly through the cav- 
ity and actually remain inside the cavity only for a short 
time. This could induce one to think that the effect 
of their presence inside the cavity might be negligible. 
We believe, however, that a comparison with a micro- 
scopic master equation describing the coupling of the en- 
tire atom-cavity system with a reservoir of electromag- 
netic modes at T temperature is highly desirable and it 
may both provide a justification of the validity of the 
phenomenological model under the typical experimental 
conditions and put into evidence the physical contexts 
where the use of such a model may be unjustified. This 
is the main motivation of the results we will describe in 
the following sections. 

We begin, in the next section, by presenting the gen- 
eral formalism, reviewed, e.g., in Ref. [l5||, to derive a 
master equation for the open quantum system of interest 
starting from the microscopic Hamiltonian of the total 
closed system (system-f environment). We will then com- 
pare the master equation we obtain with the one given 
by Eq. dSD. 



III. MICROSCOPIC DERIVATION OF THE 
QUANTUM MASTER EQUATION 

A. The general formalism 

We assume that the open quantum system of interest, 
e.g., the atom-cavity system, is part of a larger system 
whose dynamics is unitary and governed by the Hamil- 
tonian H . The external environment is that part of the 
total closed system other than the system of interest. The 
Hamiltonian of the total closed system is given by 



H = Hs + H 



E 



Hi, 



(6) 



where Hs and He, are the system and environment 
Hamiltonians, respectively, and -ffjnt is the system- 
environment interaction Hamiltonian which is taken to 
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be of the form 



Hint = A®E, 



(7) 



with A ^ and E ~ E'^ Hermitian operators acting 
on the system and on the environment Hilbert spaces, 
respectively. 

We expand the interaction Hamiltonian Hint by means 
of the relations 



A = ^AH, 



(8) 



The RWA, instead, is valid as long as the relaxation time 
of the system is much longer than the typical timescale of 
the free evolution of the quantum system, i.e., when the 
maximum of the rates 7(w) is much smaller than the min- 
imum difference between the Bohr frequencies relative to 
Hs El: 



7i„ax < Aw 

11 



(12) 



In the rest of this section we will apply this general for- 
malism to the description of cavity losses in the Jaynes- 
Cummings model. 



where n(e) is the projector onto the eigenspace corre- 
sponding to the eigenvalue e of the operator Hs and the 
sum is taken over all the Bohr frequencies relative to Hs- 
Following the standard procedure, i.e., writing down 
the Liouville-von Neumann equation for the total density 
operator in the interaction picture with respect to Hs + 
He, performing the Born-Markov and the rotating wave 
approximations, tracing out the environmental degrees of 
freedom and then going back to the Schrodinger picture, 
one obtains the following master equation for the reduced 
density operator p of the system 



quat: 



+ J2 li^) U(c.)pWAt(c.) - i {AHu;)A{c.), Pit)} 



UJ>0 



UJ>0 



A\co)p{t)A{uj)~-{A{Lo)AHu;\p{t)} 



where the relation A(— w) = A^{u!) has been used and 
where we have neglected the renormalization term [Tsj . 
The coefficients 7(w) are given by the Fourier transform 
of the correlation functions of the environment: 



+ 00 



j{u;) = / dTc'^^ {E\t)E{0)) 



(10) 



where the environment operators are in the interaction 
picture. 

We stress once more that the formalism presented 
above is valid as long as we can perform three approxima- 
tions: weak coupling or Born approximation, Markovian 
approximation and rotating wave approximation. It is 
worth recalling that the Markovian approximation can 
be seen as a coarse-graining in time, and therefore holds 
as long as the correlation time of the reservoir r is much 
smaller than the characteristic time scale of the system 
dynamics t 



B. Application of the general formalism to the 
Jaynes-Cummings model 

We model the environment as a collection of quantum 
harmonic oscillators in thermal equilibrium at T temper- 
ature and we assume that the interaction Hamiltonian 
is linear in both the electromagnetic field of the cavity 
mode and the position operators of the harmonic oscilla- 
tors, i.e., 

Hs = Hjc, HE = J2iUkb^b, 

k 

Hint = {a + a^)J29k{bk + bl}j , (13) 
fc 

with ujk the frequencies of the environment oscillators, bj, 
(bk) the creation (annihilation) operator of quanta in the 
k-th environmental mode, and gk the coupling constants. 

For this system, the operators A(w), defined in Eq. ([8]), 
are given by 

A{En>j — En^„i) = \E]\f^„i) (i?7V, m|(a + a^)|i?Ar'j) {En'jI 

= yN,N'-i(/NTl + lmVN^\EN,m) {En+i,i\, (14) 



for > 1 and 



A{E,,±-Eo)^ ^\Eo) {E,,±\ 



(15) 



t > T. 



(11) 



for = 1. In Eq. we indicate the states \En,±) 
by \En^±i) and the energy eigenvalues En^± by Eiq^±i. 
Accordingly I and m take the values ±1. 

Having in mind Eq. ([9]) we can write the Markovian 
RWA master equation for the JC model interacting with 
a thermal bath at T temperature as follows 
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En.to) 



l,m=±l N=l 
1 



{\En+i, i) (-Eat+i, i I , p} 



E 



Lm=±l N=l 



^ 2{\-^N,m) {EN,m\ , P} 



(16) 



The Kubo-Martin-Schwinger condition [if 



7(^^), 



ji-uj) = cxp - 



(17) 



ensures that the stationary state reached at time t — +oo 
is the thermal state UM 



exp 



Pth 



( Hjc 



(18) 



as expected from statistical mechanical considerations 
and as can be easily justified by the detailed balance 
principle and by means of Eq. (|17p . 

It is worth underlining a first difference between our 
microscopic master equation, given by Eq. (jl6p . and the 
phcnomenological one given by Eq. ([5]). While the ther- 
mal stationary state predicted by the phcnomenological 
model is given by 



ph 
Pth = 



exp 



(- 



Tr{exp(-^^^)}' 



(19) 



our microscopic approach predicts that the stationary 
state is the one given by Eq. (fTSl) . which differs from 
Eq. (|19p for the presence of the interaction energy term 
in Hjc [See Eq. 

We conclude this section noting the limit of validity of 
our microscopic master equation. From Eqs. ([3]) and (fT2|) 
we deduce that the RWA we have performed is valid as 
long as the smallest difference between Bohr frequencies 
relative to Hjc is much larger than the highest decay 
rate of the system, i.e. 



2n > 7n 



(20) 



since the typical evolution timescale of the system is given 
by the inverse of the Rabi frequency 20. 



IV. COMPARISON BETWEEN THE TWO 
MASTER EQUATIONS 

A. The phenomenological master equation in the 
dressed-state approximation and its relation with 
the experiments 

In the first part of this section we examine an im- 
portant feature of the phenomenological master equation 
given by Eq. ([1]) which explains its success in describing 
accurately most of the cavity QED experiments. 

In [m, [13, it been shown that, in a regime very 
close to the one in which our microscopic master equation 
is valid, i.e. for 7 ^ fi, one can approximate the phe- 
nomenological master equation by means of the so-called 
dressed-state approximation. In a later paper [8] the va- 
lidity of this approximation was carefully analyzed and it 
was discovered that it actually requires the stronger con- 
dition 7 <C ft/{2N^^^) in order to be valid. The dressed- 
state approximation amounts at neglecting, in the in- 
teraction picture with respect to the Hamiltonian Hjc, 
all the time-dependent terms oscillating at frequencies 
which are multiple of fi, under the hypothesis that this 
frequency is much larger than the decay rate 7. 

In the Schrodinger picture, the set of coupled differen- 
tial equations for the matrix elements, in the dressed- 
state approximation, {En^+\ p\En^^), (En^^\ p\En,^) 
and (-BAr.il p IEn^z^z) coincide with the corresponding set 
of equations obtained from our master equation (|16p at 
zero temperature, when the spectrum of the environment 
is flat, i.e. in the case of white noise. In other words, our 
microscopic approach justifies the validity of the dressed- 
state approximation in terms of a microscopic system- 
reservoir interaction model, and explains the success of 
the phenomenological master equation in fitting the ex- 
perimental data in the strong coupling regime, i.e., when 
the relaxation time is much longer than the frequency of 



5 



the Rabi oscillations Q. It is worth stressing, however, 
that if the spectrum of the environment is not flat the 
predictions of the two master equations differ also in the 
limit of weak damping since, in this case, the differential 
equations for the relevant matrix elements obtained from 
our microscopic master equation do not coincide with the 
differential equations obtained from the phenomenolog- 
ical master equation after performing the dressed-state 
approximation pTl [T3 | . 

In the next subsection we present a comparison be- 
tween the predictions of the phenomenological and mi- 
croscopic master equations at zero temperature when 
the system has one initial excitation. We will concen- 
trate on a finite (three) dimensional subspace using the 
phenomenological model without the dressed-state ap- 
proximation, because the restriction of Eq. ([T]) to such a 
subspace gives rise to an exactly solvable dynamical case. 

Since our microscopic model gives the same predictions 
of the dressed-state approximation for the phenomeno- 
logical model, the examples we are going to consider in 
the next subsection will give us also insight in the limits 
of validity of the dressed-state approximation itself. 

As we will see in the following, the discrepancy between 
the phenomenological master equation (without dressed- 
state approximation) and the microscopic master equa- 
tion is indeed appreciable already at the first order in 
7/f2. Consequently, for the value of the parameters con- 
sidered in our examples, the applicability of the dressed 
state approximation to the phenomenological model ap- 
pears to be questionable. 
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FIG. 1: Populations Po,g{t) and Po^g{t) vs r = 2Qt, when the 
system starts from the state |0, e), with 7/2f2 = 0.1. The solid 
line refers to the predictions given by the phenomenological 
master equation, while the dashed line refers to the predic- 
tions given by the microscopic one. In the inset it is shown the 
long time behavior of the same quantities, for < t < 100. 



tion: 



P = -i[Hjc,p] 



+ j{u;o + n) [-\Eo){Ei.+ \p\Ei^+){Eo\ 



+ j{u;o-n)[-\Eo){Ei,-\p\Ei,_){Eo\ 



(21) 



B. Dynamics at T = with one initial excitation 



1. Decay of the Rabi oscillations with the atom initially 
excited 



We assume the initial state of the system is |0, e). In 
absence of cavity losses one would observe a continu- 
ous exchange of energy between the atom and the cavity 
mode, namely the Rabi oscillations. The interaction be- 
tween the cavity and the environment causes the loss of 
energy from the atom-cavity system to the external en- 
vironment. Since in the system there is only one initial 
excitation and the environment is at zero temperature, 
the number of excitations cannot increase in time and 
Eq. (|16p reduces to the following simplified master equa- 



obtained from Eq. (fTB]) neglecting all the terms which do 
not contribute to the evolution of the system. 

In the Appendix we provide the solution of both Eq. ([l|) 
and Eq. (PTjl based on the method of the damping ba- 
sis [1^. In Fig. 1 we plot the time evolution of the 
population of the state \0,g) predicted by both the phe- 
nomenological and the microscopic master equations, 
with 7 [loq — rj) = 7 (wq + ^) = "f = ri/5. This condition 
legitimates the use of our RWA master equation, while it 
does not assure the validity of the dressed-state approx- 
imation on the phenomenological model, in accordance 
to |8i]. The populations are given by 

Po,g{t)^{0,9\p{t)\0,g)^l-e~i\ (22) 

for the predictions of Eq. (I^T]) (microscopic model), as 
one can see from Eq. (|35p in the Appendix putting 7^ = 
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7h = 7: and by 



Pt,{t)^{0,9\p{t)\0,g) 



i6r22 - ^2 



e 2- 



2(16f}2-^2) 



7V7^ 



-7-\/T2-16n2 ^ 



2(16172-72) 



(23) 



for the predictions of Eq. ^ (phenomenological master 
equation) . 

In both cases, the ground state population of the atom- 
cavity system increases in time with the same exponen- 
tial rate, due to the cavity losses. There is, however, 
an important difference in the behavior predicted by the 
two equations. Indeed, while our microscopic master 
equation predicts a purely exponential increase, the phe- 
nomenological master equation predicts the presence of 
oscillations at the Rabi frequency superimposed to the 
exponential increase. As anticipated, the amplitude of 
these oscillations is of the order of j/fl. This difference 
in the ground state population dynamics reflects different 
physical mechanisms in the dissipation process. Accord- 
ing to the phenomenological master equation, given by 
Eq. IT]), only the cavity can directly lose excitations, as 
one can see from the form of the dissipator. The rate of 
loss is indeed proportional to the population of the state 
\l,g) only, while the state |0,e) does not directly decay. 
In this sense the oscillating behavior of the population of 
the ground state is a signature of the Rabi oscillations 
which induce the system decay via the coupling of the 
state |0,e) to the state |1,5). 

The microscopic master equation given by Eq. (|16p . 
with 7 {loq — 17) = 7 (ujq + fl), predicts instead that both 
the states \1, g) and |0, e) decay with the same rate, since 
they are both superpositions of the states and 
|£'i._). For this reason there is no oscillating behavior 
in the time evolution of the population of \0,g). 

The difference between the predictions of the phe- 
nomenological master equation and those of the micro- 
scopic master equation can be revealed by performing 
joint measurements on both the atom and the cavity, in 
order to get the population of the state \0,g) of the com- 
posite system. 

We can observe a difference in the system dynamics 
also measuring the population of the atomic ground state, 
given by Pg{t) - {0, g\ p{t) \0, g) + {1, g\ p{t) \l, g) . This 
type of measurement can be performed with standard 
techniques [5] . With the help of Eq. (PT|) one finds [See 
Eqs. ^ and Eq. ^ in the Appendix] 



1 



1 r 



,(24) 



2 4 

while using Eq. (P) one finds, from Eq. and Eq. (I37| 
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FIG. 2: Populations Pg{t) and P^'^{t) vs r = 2Qt for the 
initial state |0, e), with = 0.1. The solid line refers to the 

predictions given by the phenomenological master equation, 
while the dashed line refers to the predictions given by the 
microscopic one. In the large panel, where one can see the 
long time behavior, the two lines are indistinguishable, while 
in the inset we have emphasized the dynamics for 21 < r < 26, 
to put into evidence the presence of the frequency shift. 



in the Appendix, 



16f72 - 72 



re 2' 



272 + 27V/72 - mn^ - i6f72 -^+Vt^-i6q^ 



4(16^22-^2) 
272 - 27^72 - 16^22 - 1602 -,^V-,2-ien2 ^ 

4(16172-72) ^ 



.(25) 



In this case, however, the behavior predicted by the 
phenomenological master equation is very similar to the 
one predicted by the microscopic one, as one can see from 
Fig. 2. The only difference, indeed, is the presence of 
a frequency shift in the Rabi oscillations, as predicted 
by the phenomenological model. However, since such a 
shift is of the order of (j/il)'^ [10], the differences in the 
predicted behavior may be difficult to detect. 



2. Decay of a Bell state: exponential vs. oscillatory 
behavior. 

We now consider the case in which the atom-cavity sys- 
tem is initially prepared in a Bell state, such as the state 
as given by Eq. Contrarily to the case consid- 
ered in the previous subsection, this time the measure- 
ment of the atomic ground state population Pg{t) would 
allow to bring to light the differences in the predictions of 
the phenomenological and microscopic master equations. 
For the case of the initial state \Ei+) the ground state 
population is given by 



Pg it) = 1 



(26) 
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FIG. 3; Populations Pg{t) and Pf^(t) vs r = 2Q.t, when the 
system starts from the state with "f/2Q, = 0.1. The 

sohd hne refers to the predictions given by the phenomenolog- 
ical master equation, while the dashed line refers to the pre- 
dictions given by the microscopic one. In the inset one can see 
the long time behavior of the same quantities for < r < 100. 



for the microscopic model [See Eq. (fT6|) ] and 



Pf\t) = 1 



re 2' 



+ 



4(16172-72) 

72 - 7-^/72 - 16^2 _^_v-^2_ieQ2 ^ 



4(16^2-72) 



(27) 



for the phenomenological model [See Eq. ([T])]. Details 
on the derivation of the above formulas are given in the 
Appendix. 

As we see from Fig. 3, the atomic ground state popu- 
lation, as predicted by the microscopic master equation 
given by Eq. (|16[) , exhibits a purely exponential behavior 
while the phenomenological master equation, given by 
Eq. ([T]), predicts the presence of oscillations at the Rabi 
frequency with amplitude of the order of 7/17. In other 
words, in the microscopic approach, the atom-cavity sys- 
tem as a whole is subjected to cavity losses and therefore 
the atom dressed by the cavity mode can directly decay, 
although it is only the cavity which is directly coupled to 
the environment. 

We stress that, although the generation of a Bell state 
may be more complicated than the generation of a Pock 
state, the set of measurements necessary to distinguish 
between the predictions of the two master equations 
when the initial state is is simpler than the set 

of joint measurements needed in the case of the initial 
state |0, e). In fact in the case of the Bell state we just 
need a way to experimentally distinguish between the 
atomic excited and ground state, for example by means 
of a static electric field which ionizes the excited state 
only, as shown in Q. 



DISCUSSION AND CONCLUSIVE 
REMARKS 



In the previous section we have seen that our deriva- 
tion of the dissipative dynamics of the JC model with 
losses, taking into account the presence of the atom in- 
side the cavity, allows to give a microscopic description of 
the interaction between the atom-cavity open quantum 
system and the external electromagnetic field. 

It appears evident from our results that, although the 
source of the losses is the escape of photons through the 
cavity mirrors, it is the atom-cavity system as a whole, 
i.e. the atom dressed by the cavity mode, which comes 
into play in the dissipative dynamics. Stated another 
way, while the jump operators appearing in the phe- 
nomenological master equations, given by Eq. ([1]), are 
the annihilation and creation operators of photons in the 
cavity (a, and a^), those appearing in the microscopic 
master equation [See, e.g., Eq. (pij) ] describe transitions 
between the dressed states of the atom-cavity system. 

A noteworthy point emerging from our analysis is that 
the deviations of the exactly solved phenomenological 
model from our microscopic model (coinciding with the 
one obtained from the phenomenological model with the 
dressed-state approximation) are of the first order in 7/17. 
For this reason it should be possible to check experimen- 
tally the predictions of the microscopic model with the 
set up currently used in cavity QED experiments. 

It is worth stressing that in Eq. ([T]) only one cavity loss 
rate 7 = ^{ujf)), with ^{lo) given by Eq. (fTU]). appears. 
On the contrary, in Eq. ([2T|) . the loss rates correspond- 
ing to jump operators connecting the ground state with 
the states and [i^i.-) are different, and more pre- 

cisely they are given by 7(^0 + 17) and 7(a'o — 17), respec- 
tively. This clearly indicates that, when the spectrum 
of the reservoir is not flat, the phenomenological mas- 
ter equation does not provide a description which can be 
justified in terms of a microscopic system-reservoir inter- 
action model as the one we have considered. 

In the light of the considerations made above we ex- 
pect that, already in the weak coupling regime, differ- 
ences between our approach and the phenomenological 
one should become evident in those physical contexts 
where the reservoir is structured, e.g. in photonic band 
gap materials and hybrid solid-state cavity QED sys- 
tems Slillii. 
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APPENDIX 



In this appendix we recall the method of the damping 
basis introduced in Ref. to solve master equations, 
and we apply it to Eq. (PT|) . 

Given a master equation of the form 



(28) 



where £ is a time-independent linear superoperator act- 
ing on /9, one considers the following eigenvalue problem 



(29) 



where p\ is a right eigenoperator of the superoperator C 
with eigenvalue A. When the set of right eigenoperators 
{p\\ is a basis for the space of linear operators acting on 
the Hilbert space of the system, as in the case analyzed 
in this paper, any density operator can be expanded with 
respect to this set. 

It is easy to show that if the system starts from the 
initial condition 



P(0) 



X 



(30) 



then the time evolution of the system is then given by 



p{t) 



Y^cxe^'px 



(31) 



The coefficients {cx} of the decomposition in Eq. ([30] 
are given by 



cx = Tr{pxpm. 



(32) 



where px is the solution of the left eigenvalue problem 



pxC = Xpx, 



(33) 



where the set of left eigenvalues {A} is the same as in the 
right eigenvalue problem 

When the dimension of the Hilbert space one consid- 
ers is finite and equal to n, the superoperator can be 
represented by a non-hermitian x matrix. The 
right eigenoperators, if they exist, are represented by n^- 
component column vectors, while the left eigenoperators 
are represented by n^-component row vectors. 

Applying this method to Eq. pT|) . we have to look for 
the nine right eigenoperators of its superoperator. For 
simplicity we call 7a = 7 (wq — ri) and — 1 (wq -|- fi). 

Six of the nine right eigenoperators are given by 
the coherences \Eo) {Ei^^\, \Eq) {Ei^+\, \Ei^^) {Ei^+\ 
and their hermitian conjugates, with eigenvalues re- 
spectively equal to i{u)a — Vt) — 7a/2, i(ijjQ + Vt) — 
76/2, i(2ri) — (7a+7b)/4 and their complex con- 
jugates. The other three right eigenoperators 
are given by \Eo) {E^l {\Ei,_) {Ei^_\ - \Ea) {Eq\), 
(|£^i^+) (-Ei.+ I — l-Eo) {Eo\) with eigenvalues respectively 
equal to 0, — 7a/2 and —76/2. 



By expanding the initial state with respect to these 
nine eigenoperators, one can then compute the evolution 
of the system under the initial conditions given in part 
B of Sec. IV. 

When the initial state is p(0) = |0, e) (0, e| one obtains 
the following density operator at time t 



1 



1 



|£^o) (Sol 



+ ie-** H- ie-** \E,^+) {E,^+\ 

_ le-^* (e^^^* | + h.c.) , (34) 

which allows us to express the population of the state 
|0, as follows 



Po,<,(i) = (0,.g|p(t)|0,.g) = l 



1 —22.+ 1 „2Ly 



2 2 

The population of the state 1 1 , 5) is given by 



(35) 



Pi,3(t) = (l,5|p(i)|l,5)-7 



e 2 



2if2- 



1a+1b 



)*_eH"-^)*' 



(36) 



Using the two equations above one can calculate Pg [t) — 
{g\p{t) \g) = {Q,g\p{t)%g) + {l,g\p{t)\l,g) given in Eq. 
dill). 

Equation ([35]) and Eq. (f36|) are to be compared with 
the corresponding quantities computed using the phe- 
nomenological master equation, given by Eq. ([T]), and 
considering the same initial condition. Following the 
same lines one derives Eq. (1^^ for the population of the 
ground state |0,(?) and 



PtAt) = {l,g\p{t)\l,g) = 



16f72 - 72 



e 2 



4(16f72-72) 
161^2 

'4(16172-72) 



(37) 



for the population of state |1,5). 

In the same way one can derive the time evolution 
predicted by the microscopic master equation ([2T|) when 
the initial condition is p(0) = (-^i.+h 



p{t) =(i-e-**)|£;o)(£;o| 



+c-^^\E,,+){E,M, (38) 
which gives the following expressions for the populations 
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of states \0,g) and \l,g) for the population of \0,g), and 
Po,9W = (0,5|pW|0,<?>-l-e-^*, 

Pi^,W = (l,5|p(i)|l,<?)-^c-*'. (39) 

^f.^W = (1,.9IpW|1,.9> = 77^7^— le-i* 

Once more, we must compare these expressions with " ' 

the corresponding ones obtained by means of the phe- 2 /~2 ^5^2 ^ 2^2 

nomenological master equation given by Eq. Using _i 2X2 e * 

this equation one obtains ^ (1651^ — 7^) 



1602 ^ 7^+7^7^ - 16f^^2 -y-V-,-^-i6u^ 



7^ -^ + V72_i6q2 ^ 

-I ■ e 2 

2(16f72-72)'^ 



7^ -2-V2^li^+ 
+ 2{16n^-j-^f ' ' ^"^^^ for the population of I 
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